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Abstract 



We calculate the critical value of the hopping parameter, k c , in 0(a) 
improved Lattice QCD, to two loops in perturbation theory. We employ the 
Sheikholeslami-Wohlert (clover) improved action for Wilson fermions. 

The quantity which we study is a typical case of a vacuum expectation 
value resulting in an additive renormalization; as such, it is characterized by 
a power (linear) divergence in the lattice spacing, and its calculation lies at 
the limits of applicability of perturbation theory. 

The dependence of our results on the number of colors N, the number of 
fermionic flavors Nf, and the clover parameter cs\y> is shown explicitly. We 
compare our results to non perturbative evaluations of k c coming from Monte 
Carlo simulations. 
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I. INTRODUCTION 



In this paper we calculate the critical value of the hopping parameter k c in Lattice 
QCD, to two loops in perturbation theory. We employ the 0(a) improved Sheikholeslami- 
Wohlert (clover) action for Wilson fermions; this action is widely used nowadays in Monte 
Carlo simulations, as a means of reducing finite lattice spacing effects, leading to a faster 
approach to the continuum. 

The Wilson fermionic action is a standard implementation of fermions on the lattice. It 
circumvents the notorious doubling problem by means of a higher derivative term, which 
removes unphysical propagator poles and has a vanishing classical continuum limit; at the 
same time, the action is strictly local, which is very advantageous for numerical simulation. 
The price one pays for strict locality and absence of doublers is, of course, well known: The 
higher derivative term breaks chiral invariance explicitly. Thus, merely setting the bare 
fermionic mass to zero is not sufficient to ensure chiral symmetry in the quantum contin- 
uum limit; quantum corrections introduce an additive renormalization to the fermionic mass, 
which must then be fine tuned to have a vanishing renormalized value. Consequently, the 
hopping parameter k, which is very simply related to the fermion mass, must be appropri- 
ately shifted from its naive value, in order to recover chiral invariance. 

By dimensional power counting, the additive mass renormalization is seen to be linearly 
divergent with the lattice spacing. This adverse feature of Wilson fermions, typical of vacuum 
expectation values of local objects, poses an additional problem to a perturbative treatment, 
aside from the usual issues related to lack of Borel summability. Indeed, our calculation 
serves as a check on the limits of applicability of perturbation theory, by comparison with 
non perturbative results coming from Monte Carlo simulations. 

In the present work we will follow the procedure and notation of Ref. ||, in which k c 
was computed using the Wilson fermionic action without 0(a) improvement. The results of 
Ref. were recently confirmed in Ref. ||, in which a coordinate space method was used to 
achieve even greater accuracy. 

The critical fermionic mass and hopping parameter will now depend not only on the 
number of colors N and of fermionic flavours Nf, but also on the free parameter csw which 
appears in the clover action (see next Section); we will keep this dependence explicit in our 
results. 

In Sec. |T|we define the quantities which we set out to compute, and describe our calcula- 
tion. In Sec. |H we present our results and compare with Monte Carlo evaluations. Finally, 
in Sec. [TV] we obtain improved estimates coming from a tadpole resummation. 



Our starting point is the Wilson formulation of the QCD action on the lattice, with the 
addition of the clover (SW) |L| fermion term. Its action reads, in standard notation: 



II. FORMULATION OF THE PROBLEM 
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- ^ E E l^f( x ) ( r ~ 7 M ) U h{ x )^A x + A) + ^f(x + p.) (r + 7m ) U^x)^ f {x) 

Z f x,n 

1 _ ^ 

+ 7CswE E ^f{x)a^F^{x)-lj} f {x), (1) 
f 

where: = ^ (Q^ - Q v/X ), = U^ v + C/^ _ M + C/^, _„ + [/_„, ^ (2) 

o 

Here U^ u (x) is the usual product of link variables U^(x) along the perimeter of a plaquette 
in the \i-v directions, originating at x; go denotes the bare coupling constant; r is the Wilson 
parameter; / is a flavor index; = (i/2)[ / y fJi , j v ]. Powers of the lattice spacing a have been 
omitted and may be directly reinserted by dimensional counting. 

We use the standard covariant gauge-fixing term; in terms of the vector field Q^x) 
[U^(x) = expiigoQ^x))], it reads: 

Sgf = Ao E E Tr &-Q»(x)A;Q u (x), A-Q v (x) =Q v (x-p,)- Q u (x). (3) 

fl,U X 

Having to compute a gauge invariant quantity, we chose to work in the Feynman gauge, 
A = 1. Covariant gauge fixing produces the following action for the ghost fields u and To 

[Q^x),[Q^x),A+u(x)]]+---), A+u(x)=u(x + p,)-u(x). (4) 



x [J, 
12^0 



Finally the change of integration variables from links to vector fields yields a jacobian that 
can be rewritten as the usual measure term S m in the action: 

S m = ^Ngl EE Tr WW + • • • (5) 

In S g h and S m we have written out only terms relevant to our computation. The full action 
is: S = Sl + S g f + S g h + S m . 

The bare fermionic mass m# must be set to zero for chiral invariance in the classical 
continuum limit. The value of the parameter Csw can be chosen arbitrarily; it is normally 
tuned in a way as to minimize 0(a) effects. Terms proportional to r in the action, as well 
as the clover terms, break chiral invariance. They vanish in the classical continuum limit; at 
the quantum level, they induce nonvanishing, flavor-independent corrections to the fermion 
masses. Numerical simulation algorithms usually employ the hopping parameter, 

(6) 



2 rriB a + 8 r 

as an adjustable quantity. Its critical value, at which chiral symmetry is restored, is thus 
l/8r classically, but gets shifted by quantum effects. 

The renormalized mass can be calculated in textbook fashion from the fermion self- 
energy. Denoting by m B , g ) the truncated, one particle irreducible fermionic two- 
point function, we have for the fermionic propagator: 
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S(p) = [if + m(p) - S L (p, m B , g )] ' (7) 

1 2r 
where : § = 7^— sin(ap M ), m(p) = H sin 2 (ap M /2). 



M a a M 



To restore the explicit breaking of chiral invariance, we require that the renormalized 
mass vanish: 

,S- 1 (0) = m B = E L (0,m B , 5o ) (8) 

The above is a recursive equation for mg, which can be solved order by order in perturbation 
theory. 

We write the loop expansion of S L as: 

£ L (0, m B , g ) = gl £« + g 4 + • • • (9) 

Two diagrams contribute to £W, shown in Fig. I. In these diagrams, the fermion mass 
must be set to its tree level value, m B — > 0. 




1 2 

FIGURE I. One-loop diagrams contributing to Y1 L . Wavy (solid) lines represent gluons (fermions). 

The quantity S^ 2 ^ receives contributions from a total of 26 diagrams, shown in Fig. II. 
Genuine 2-loop diagrams must again be evaluated at itlb — > 0; in addition, one must include 
to this order the 1-loop diagram containing an O(g^) mass counterterm (diagram 23). 

The contribution of the i th diagram can be written in the form 

(^-D-t(cg + f + ^c«Aw (10) 

where cfj, c^], c^J are numerical constants. The dependence on csw is seen to be polynomial 
of degree 4, as can be verified by inspection of Fig. II. 

Certain sets of diagrams, corresponding to renormalization of loop propagators, must 
be evaluated together in order to obtain an infrared-convergent result: These are diagrams 
7+8+9+10+11, 12+13, 14+15+16+17+18, 19+20, 21+22+23. 
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FIGURE II. Two- loop diagrams contributing to Y1 L . Wavy (solid, dotted) lines represent gluons 
(fermions, ghosts). Crosses denote vertices stemming from the measure part of the action; a solid 
circle is a fermion mass counterterm. 



III. NUMERICAL RESULTS 

Evaluating the two diagrams of Figure I, we find for : 

E (i) = N 2 - 1 ( _0.15493339023106021 (diagram 1) 

-0.00792366847979(1) +c sw 0.04348303388205(10) 

+c sw 2 0.01809576878142(1) ) (diagram 2). 

(11) 

Here and below we set r to its usual value, r = 1 . One- and two-loop results pertaining to 
c sw — are as in Ref. @, and can be found with greater accuracy in Ref. @. 



5 



For csw 7^ 0, only one-loop results exist so far in the literature; a recent presentation for 
the case csw=l ( see Ref. @] an d earlier references therein) is in perfect agreement with our 
Eq. (0): 



£«(iV=3,csw=l) = -0.2700753495(2) [Ref. [4]] 

-0.2700753494597(5) [The present work, Eq. (|TJ)]. ^ > 

It should be clear to the reader that such a high level of precision is hardly relevant 
per se, especially given the expected deviation from nonperturbative results; nevertheless, 
it serves as a testing ground for both accuracy and efficiency of our methods, in view of the 
more demanding higher loop calculations. Further, high precision is called for in the context 
of the Schrodinger functional computation, to permit a stable extrapolation of various pa- 
rameters to infinite lattice (see Ref. @). Regarding efficiency, let us note that the numerical 
integrations leading to Eq. (O) require a mere ~ 10 min of CPU time on a typical 1GHz 
Pentium III processor. 

We now turn to the much more cumbersome evaluation of the two-loop diagrams, which 
is the crux of our present computation. As in Ref. ||, we use a Mathematica package which 
we have developed for symbolic manipulations in lattice perturbation theory (see, e.g., Ref. 

For the purposes of the present work (and a related work on the /9-function 0), we 
have augmented the package to include the vertices of the clover action. 

In Tables I, II, III and IV we present the values of the coefficients cj£], c£j, c^J, c£J, 

respectively. The C(csw°) coefficients cf] are as in Ref. [0, and have been listed in Table 
V for completeness. Diagrams giving vanishing contributions to a given power of csw have 
been omitted from the corresponding table. 

The momentum integrations leading to the values of each coefficient are performed nu- 
merically on lattices of varying size L < 32, and then extrapolated to infinite lattice size 
using a broad spectrum of functional forms of the type: J2i,j % (In L) 5 / U. The systematic 
error resulting from the extrapolations has been estimated rather conservatively using the 
procedure of Ref. |5|], and has been included in the tables. 

One important consistency check can be performed on those diagrams which are sepa- 
rately IR divergent; taken together in groups, as listed at the end of Section II, they give 
finite and very stable extrapolations for the coefficients of each power of csw- Several other 
consistency checks stem from exact relations among various coefficients; to name a few: 

c« = (1/2 -&?))/ 4 



c 



(0) _^ 1,(0) 

2,7 — °1 

Jl) _ 1,(1)/. 

c 2,14 — u 2 I ' 

(2) _ _l(2) /; 

c 2,14 — u 2 I ' 



6i Uj 3/32 (13) 



Here, are the coefficients of the iih 1-loop diagram, multiplying csw j , as displayed in 



Eq. (p4|). Comparing with our numerical values of Tables I- IV, we find agreement well within 
the error bars. 

Leaving the choice of values for N, Nf and csw unspecified, our result takes the form: 
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£(2) = (N 2 -l)[ 

+ 
+ 
+ 
+ 



-0.017537(3) +1/N 2 0.016567(2) +N f /N 0.00118618(8)) 

0.002601(2) -1/N 2 0.0005597(7) -N f /N 0.0005459(2)) c sw 

-0.0001556(3) +1/N 2 0.0026226(2) +N f /N 0.0013652(1)) c sw 2 

-0.00016315(6) +1/N 2 0.00015803(6) -N f /N 0.00069225(3)) c sw 3 

-0.000017219(2) +1/N 2 0.000042829(3) -N f /N 0.000198100(7)) c sw 4 ] 

(14) 



To make more direct contact with non-perturbative results, we evaluate £( 2 ) at iV = 3 
and Nf = 0, 2, obtaining: 

£( 2 )(jV = 3,N f = 0) = -0.12557(3) +0.02031(2) c sw +0.001087(3) c sw 2 

-0.001165(1) c sw 3 -0.00009968(2) c sw 4 

£( 2 )(jV = 3, N f = 2) = -0.11924(3) +0.01740(2) c sw +0.008368(3) c sw 2 ' 

-0.004857(1) c sw 3 -0.0011562(1) c sw 4 

Eqs. ( [II] , [T5D lead immediately to the 1- and 2-loop results for the critical mass: = 
^pX 1 - 1 ), m^P= + c?qS^ 2 ' 1 , and the corresponding hopping parameter k c = 1/(2 m c a + 

8r). 

A number of non-perturbative determinations of n c exist in the literature for particular 
values of /3 = IN j g\ and csw 7^ 0, see e.g. Refs. |7||| (quenched case) and Refs. PJlQfl 
(unquenched, Nf = 2). We present these in Table VI, together with the 1- and 2-loop 
results (^4 , Also included in the Table are the improved results obtained with the 

method described in the following Section. 



IV. IMPROVED PERTURBATION THEORY 

In order to obtain improved estimates from lattice perturbation theory, one may perform 
a resummation to all orders of the so-called "cactus" diagrams [pT|-|l3|. Briefly stated, these 
are gauge-invariant tadpole diagrams which become disconnected if any one of their vertices 
is removed. The original motivation of this procedure is the well known observation of 
"tadpole dominance" in lattice perturbation theory. In the following we adapt the calculation 
of Ref. P| to the clover action. We refer to Ref. [JTTJ] for definitions and analytical results. 

Since the contribution of standard tadpole diagrams is not gauge invariant, the class 
of gauge invariant diagrams we are considering needs further specification. By the Baker- 
Campbell-Hausdorff (BCH) formula, the product of link variables along the perimeter of a 
plaquette can be written as 

= exp {ig (A Xj)X + A x+ ^ v - A x+V ^ - A Xjl/ ) + 0(go)} 

= exp {igoFi%, + ig 2 F^ v + O(g 4 )} (16) 

The diagrams that we propose to resum to all orders are the cactus diagrams made of 
vertices containing F^ u . Terms of this type come from the pure gluon part of the lattice 
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action. These diagrams dress the transverse gluon propagator Pa leading to an improved 
propagator pjj , which is a multiple of the bare transverse one: 

PT = r , (17) 

where the factor w(go) will depend on g and N, but not on the momentum. The function 



w(go) can be extracted by an appropriate algebraic equation that has been derived in Ref. [3 I 
and that can be easily solved numerically; for SU(3), w(go) satisfies: 



ue~ u/3 



M 2/ 3 _ Au + 8 ] = 2g 2 , u(g ) = 9 % ■ (18) 



The vertices coming from the gluon part of the action, Eq. ([I]), get also dressed using a 
procedure similar to the one leading to Eq. dl7|) [fllf| . Vertices coming from the Wilson 
part of the fermionic action stay unchanged, since their definition contains no plaquettes on 
which to apply the linear BCH formula; the 3- and 4-point vertices of the clover action, on 
the other hand, acquire simply a factor of (1 — w(go)) [p^j. 

One can apply the resummation of cactus diagrams to the calculation of additive and 
multiplicative renormalizations of lattice operators. Applied to a number of cases of inter- 



est |pl] , |12| , this procedure yields remarkable improvements when compared with the avail- 
able nonperturbative estimates. As regards numerical comparison with other improvement 
schemes (tadpole improvement, boosted perturbation theory, etc.) |14|,|15l , cactus resumma- 
tion fares equally well on all the cases studied JTJ| . 



One advantageous feature of cactus resummation, in comparison to other schemes of 
improved perturbation theory, is the possibility of systematically incorporating higher loop 
diagrams. The present calculation exemplifies this feature, as we will now show. 

Dressing the 1-loop results is quite straightforward: the fermionic propagator stays un- 
changed, the gluon propagator gets multiplied by 1/(1 — w(go) and the dressing of the 
fermionic vertices amounts to a rescaling: csw — > c sw (1 — w (9o))- The resulting values, 
^dressed) are shown in Table [V]]. It is worth noting that these values already fare better 



than the much more laborious undressed 2-loop results. 

We now turn to dressing the 2-loop results. Here, one must take care to avoid double 
counting: A part of diagrams 4, 7 and 14 has already been included in dressing the 1-loop 
result, and must be explicitly subtracted from E^ 2 ) before dressing. Fortunately, this part 
(we shall denote it by E^y is easy to identify, as it necessarily includes all of the 1/iV 2 
part in diagrams 7, 14, and the 1/iV 2 part of diagram 4 involving a clover 5-point vertex. 
A simple exercise in contraction of SU(N) generators shows that is proportional to 
(2iV 2 - 3)(iV 2 - l)/(3iV 2 ). There follows without difficulty that: 

S^ b = -(2iV 2 -3)(iV 2 -l)/(3iV 2 )- 

[t>2 Csw/o + C 2>4 CSW + C 2i7 + C 214 + C 2 U C S W + C 2 ,14 C SW J {^) 

(cf. Eq. (0)). 

A potential complication is presented by gluon vertices. While the 3-gluon vertex dresses 
by a mere factor of (1 — w(g )), the dressed 4-gluon vertex contains a term which is not 
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simply a multiple of its bare counterpart (see Appendix C of Ref. ||11||). It is easy to check, 
however, that this term must simply be dropped, being precisely the one which has already 
been taken into account in dressing the 1-loop result; the remainder dresses in the same way 
as the 3-gluon vertex. The very same situation prevails with the dressed 5-point vertex of 
the clover action, as in diagram 4. 

In conclusion, cactus resummation applied to the 2-loop quantity Y,^ leads to the fol- 
lowing rather simple recipe: 



rn 



(2) 

c, dressed 



E (D 



1 - w(g 



( E (2) 



y (2) s gt 

^sub/ fl „../'„ \"|2 



[i -w(go) 



(20) 



csw^cswCi-^Cflo)) 



(For the values of (3 in Table 0, (3 = 5.20, 5.26, 5.29, 5.7, 6.0, 6.2, 12.0, 24.0, we obtain from 
Eq. QT3): l-w(g ) = 0.697146, 0.701957, 0.704298, 0.732579, 0.749775, 0.759969, 0.887765, 
0.946087, respectively.) 

Our results for /^dressed) as obtained from Eq. (|20"D, are listed in Table [VT|. Comparing 
with the Monte Carlo estimates, dressed results show a definite improvement over non- 
dressed values. It is interesting to note that 1-loop dressed results already provide most of 
the improvement, except at very large /3-values. At the same time, a sizeable discrepancy 
still remains, as was expected from start; multiplicative renormalizations, calculated to the 
same order, are expected to be much closer to their exact values. A first case study of this 
kind, regarding the /3-function with clover improvement, has now been completed and is 
presented in Ref. ||. 
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TABLE I. Coefficients c\ (, & 2 



(i) „M Jl) 

i i °3,i m 



„(L 
-l.i 



2,1 



-3,i 



4 

12+13 
14+15+16+17+18 
19+20 
21+22+23 
24 
25 
26 
27 
28 



-0.00697298969(1) 


0.005587(1) 


0.0015499(6) 
-0.00022742(6) 
0.0014715(3) 
0.0009439(2) 


0.00024887(3) 



0.00711962270(1) 


-0.0054356(3) 


-0.0015499(6) 


-0.00002752(1) 


-0.0007525(1) 
0.000086138(5) 





-0.00005540(1) 


-0.0004905(2) 
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TABLE II. Coefficients 


(2) (2) (2) 




i 


c i,i 


C 2,i 


J 2 ) 

c 3,i 


4 


-0.00486917062(1) 


0.00452394220(1) 





6 


0.0017538(2) 








12+13 








0.0008949(1) 


14+15+16+17+18 


0.0021977(2) 


-0.0022620(2) 





19+20 








0.0004703(1) 


21+22+23 


-0.0001864(1) 


0.0001864(1) 





24 


0.00003257(1) 








25 


-0.00022829(1) 


-0.00005915(1) 





26 


0.00060875(1) 








27 





0.00035168(2) 





28 


0.00053539(5) 


-0.00011818(1) 







TABLE III. Coefficients 


J3) _(3) (3) _ 1 

1,2 " " o,2 




i 


(3) 


(3) 


(3) 


19+20 








-0.00069225(3) 


21+22+23 


-0.00017530(6) 


0.00017530(6) 





25 


0.000022090(4) 








26 


-0.000023954(3) 








27 





-0.000017264(2) 





28 


0.000014014(3) 










TABLE IV. Coefficients 


J4) J4) (4) 1 

c l,i ' c 2,i > c 3,i ■ r — l. 




i 


c (4) 


c (4) 


c {4) 


19+20 








-0.00019810(1) 


21+22+23 


-0.000017219(2) 


0.000017219(2) 





27 





0.000025610(2) 
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TABLE V. Coefficients c$, 6%!, cgj. r = 1. 



z 


c (0) 


c (0) 


c {0) 


3 


0.002000362950707492 


-0.0030005444260612375 





4 


0.00040921361(1) 


-0.00061382041(2) 





6 


-0.0000488891(8) 


0.000097778(2) 





7+8+9+10+11 


-0.013927(3) 


0.014525(2) 





12+13 








0.00079263(8) 


14+15+16+17+18 


-0.005753(1) 


0.0058323(7) 





19+20 








0.000393556(7) 


21+22+23 


0.000096768(4) 


-0.000096768(4) 





25 


0.00007762(1) 


-0.00015524(3) 





26 


-0.00040000(5) 








27 





-0.000006522(1) 





28 


0.0000078482(5) 


-0.000015696(1) 






TABLE VI. One- and two-loop results (kc , n c ), along with their improved (dressed) 
counterparts, and nonperturbative determinations. See references, shown in square brack- 
ets, for details on the nonperturbative definition of n and on error estimates. 



N f 


P 


csw 










JXi 

c, dressed 


K {2) 
c, dressed 




Simulation 





5.70 


1.568 





1296 





1332 


0.1366 


0.1366 





1432 







6.00 


1.479 





1301 





1335 


0.1362 


0.1362 





1392 


§ 





6.00 


1.769 





12749 





13061 


0.13372 


0.13319 





13525 







6.20 


1.442 





1303 





1334 


0.1358 


0.1358 





1379 


§ 





6.20 


1.614 





12878 





13182 


0.13439 


0.13414 





13582 







12.0 


1.1637 





128766 





129622 


0.129807 


0.129845 





129909 








24.0 


1.0730 





127019 





127229 


0.127243 


0.127253 





127258 


a 


2 


5.20 


2.0171 





12515 





12987 


0.13481 


0.13342 





13663 


|1C,| 


2 


2.26 


1.9497 





12589 





13043 


0.13517 


0.13392 





13709 


MB 


2 


2.29 


1.9192 





12622 





13068 


0.13532 


0.13414 





13730 


01 
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